Vector boson production in association with KK 
modes of the ADD model to NLO in QCD at LHC 

M. C. Kumar'^Q Prakash Mathews^@ V. Ravindran''@ Satyajit Seth^Q 



" Regional Centre for Accelerator-based Particle Physics 
Harish- Chandra Research Institute, Chhatnag Road, Jhunsi, 
Allahabad 211 019, India 

^ Saha Institute of Nuclear Physics, 1/AF Bidhan Nagar, Kolkata 700 064, India 



Abstract 

Next-to-leading order QCD corrections to the associated production of vector 
boson {Z /W^) with the the Kaluza-Klein modes of the graviton in large ex- 
tra dimensional model at the LHC, are presented. We have obtained various 
kinematic distributions using a Monte Carlo code which is based on the two 
cut off phase space slicing method that handles soft and collinear singularities 
appearing at NLO level. We estimate the impact of the QCD corrections on 
various observables and hnd that they are significant. We also show the re- 
duction in factorization scale uncertainty when QCD corrections are included. 
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1 Introduction 



With the on set of the Large Hadron Colhder (LHC) era, unique opportunity to 
probe the realm of new physics in the TeV scale has begun. Models with extra 
spatial dimensions and TeV scale gravity, are proposed to address the large hierarchy 
between the electroweak and Planck scale and are expected to provide a plethora of 
new and interesting signals. 

The extra dimension model proposed by Arkani-Hamed, Dimopoulos and Dvali 
(ADD) [1], was the first extra dimension model in which the compactified dimen- 
sions could be of macroscopic size and consistent with present experiments. A viable 
mechanism to hide the extra spatial dimension, is to introduce a 3-brane with negli- 
gible tension and localise the Standard Model (SM) particles on it. Only gravity is 
allowed to propagate in the full 4 + 5 dimensional space time. As a consequence of 
these assumptions, it follows from Gauss Law that the effective Planck scale Mp in 
4-dimensions is related to the 4 -(- 5 dimensional fundamental scale Ms through the 
volume of the compactified extra dimensions [1]. The extra dimensions are compact- 
ified on a torus of common circumference Rs- The compactified extra dimensions 
are flat, of equal size and could be large. The large volume of the compactified extra 
spatial dimensions would account for the dilution of gravity in 4-dimensions and 
hence the hierarchy. Current experimental limits on deviation from inverse square 
law [2], constraint the number of possible extra spatial dimensions 5 > 2. The space 
time is factorisable and the 4-dimensional spectrum consists of the SM confined to 
4-dimensions and a tower of Kaluza-Klein (KK) modes, of the graviton propagating 
in the full 4 + 5 dimensional space time. 

The interaction of the KK modes G^i^J with the SM fields localised on the 3-brane 
is described by an effective theory given by [3l H] 



where T^'^ is the energy-momentum tensor of the SM fields on the 3-brane and 
Mp = Mpj^^i is the reduced Planck mass in 4-dimensions. The relation between 
the 4-dimensional coupling, the volume of the extra dimensions and the fundamental 
scale Mp, in 4 + 5-dimensions 



The size of the extra dimension Rs is related to the radius Rp^ Rs = 2'nRp. The 
fundamental scales in 4 + 5 dimensions, as defined in [3] Mp is related to Ms [1] 
as: 87rM^+2 = M^+^S^-i, where S5-1 = 2^^/"^ /T{5/2) is the surface area of a unit 
sphere in S dimensions. 

The zero mode corresponds to the usual 4-dimensional massless graviton and 
higher massive KK modes are labeled by n = (ni,n2, ■ ■ ■ ,ns). The masses of the 
individual KK modes are m| = •n^/R'jj and the mass gap between adjacent KK 
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modes is Am = Rj^. For not too large 5 the discrete mass spectrum could be 
replaced by a continuum, with the density of KK states 

Pinirt) = \Ss^lRWn~^ ■ (3) 

For an inclusive cross section at the collider, we have to sum over all accessible KK 
modes and hence cross section for the production of an individual KK mode of mass 
rrifi has to be convoluted with the density of states pirrifi). The discrete sum YIih 
can be replaced by J dm? p{mfi) , and hence the inclusive cross section for the tower 
of KK modes is 

d(j = Ss-i — TTTj [ dm m^^^ da^jf* , (4) 
M|+'^ J ^ ^ ' 

where dcx^'* is the cross section to produce a single KK mode. The cross section 
for an individual KK mode is suppressed by the coupling factor (2Mp)~^, the high 
multiplicity of accessible KK modes at the collider would compensate, leading to the 
exciting possibility of observing low scale quantum gravity effects at the LHC The 
additional 1/2 factor is due to the definition of the sum of polarisation of the KK 
modes [3]. The is kinematically constrained to those KK modes which satisfy 

= \^\/Rd < y/s, where ^/s is the partonic center of mass energy or as the case 
may be the available energy to produce the KK modes. 

Viable signatures of the ADD scenario at the LHC are possible by the exchange of 
virtual KK modes between the SM particles, leading to an enhanced cross section or 
by the emission of real KK modes from the SM particles, leading to a missing energy 
signal. Various such processes have been extensively studied in this model, most of 
which have been considered only up to leading order (LO) in QCD [21 IH El El E]- 
These LO approximations at the hadron colliders suffer from large factorisation and 
renormalisation scale dependence which for some processes could be as large as a 
factor of two. These issues go beyond normalisation of a cross section by a i^"-factor 
as the shapes of distributions may not be modeled correctly and in addition the 
LO cross sections are strongly dependent on the factorisation scale. It is hence 
essential to evaluate the next-to-leading order (NLO) corrections to the process 
of interest to provide quantitatively reliable predictions. NLO QCD corrections 
to extra dimension models have been studied for dilepton [8], boson pair [9l [10] 
productions, and real graviton production processes such as graviton plus jet [TT] 
and graviton plus photon [12]. Searches at the Tevatron using the single photon or jet 
with missing transverse energy have been used to put bounds on extra dimensional 
scale Md for different number of extra dimensions [13^ il4j. The same signal has 
been simulated for the LHC at the ATLAS detector [15], discovery limit and the 
methods to determination of the parameters of the extra dimensional models are 
discussed. 

In this paper we consider the graviton production in association with a vector 
boson at the hadron colliders at NLO in QCD. Z-boson process to LO had been 
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considered at LEP [16j and simulation studies for the Z + Gkk modes production 
at the LHC was studied to LO ^7] as a complement to the more conventional 
channels. 

2 Vector bosons in association with KK modes 

The associated production of Z-boson with the KK modes of the ADD model leads 
to missing energy signals at the hadron colliders. We begin by discussing the neutral 
weak gauge boson (Z) production in association with the KK modes of the ADD 
large extra dimensional model to NLO in QCD and would consider the charged weak 
gauge bosons (VT^) towards the end. 

At the hadron collider, the associated production P P Z Gkk X at LO 
proceeds via the quark, anti-quark annihilation process q q ^ Z Gkk- There are 
four diagrams that would contribute to this process, which corresponds to the KK 
modes of the graviton being emitted of a fermion leg, Z-boson or the q q Z vertex. 
The Feynman rules to evaluate the matrix elements are given in [3l H] and for the 
vector boson, unitary gauge — )■ oo) is used. Summation of the polarisation tensor 
of the KK modes is given in It can be seen that the terms proportional to the 
inverse powers of KK mode mass vanish on expressing the matrix element square 
in terms of independent variable and is an useful check. 

The NLO calculation presented here uses both analytical and Monte Carlo in- 
tegration methods and hence is flexible to incorporate the experimental cuts and 
can generate the various distributions unlike a fully analytical computation. Our 
code is based on the method of two cutoff phase space shcing |l8j to deal with 
various singularities appearing in the NLO computation of the real diagrams and 
to implement the numerical integrations over phase space. The analytical results 
are evaluated using the algebraic manipulation program FORM The real and 
virtual corrections have been evaluated in the massless quark limit. We use dimen- 
sional regularisation with space time dimensions d = 4 + e. To deal with 75 in 
(i- dimensions, we use the completely anti-commuting 75 prescription [20] . 

The order 0{as) corrections to the associated production of the Z-boson and 
KK modes of the graviton come from the following 2 — )■ 3 real diagrams (a) g g — > 
Z Gkk 9, (b) q{q) g — )■ q{q) Z Gkk- There are 14 diagrams that contribute to the 
quark antiquark annihilation process and can be classified into diagrams where the 
KK modes couples to the (i) fermion legs, (ii) Z-boson leg, (iii) gluon leg, (iv) q q Z 
vertex and {v) q q g vertex. The unitary gauge is used for the Z boson propagator 
and Feynman gauge for the gluon propagator. Note that the gauge parameters 
influence the coupling of KK modes to the SM fields For the sum over the 
polarisation vectors we have retained only the physical degrees of freedom. 

The real diagrams involving gluons and massless quarks would be singular in 
soft and coUinear regions of the 3-body phase space integration. Two small slic- 
ing parameters Ss and 6c are introduced to isolate regions of phase space that are 
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sensitive to soft and collinear singularities. Rest of the region is finite and can be 
evaluated in 4-dimensions. Phase space integrations in the mutually exclusive soft 
and collinear regions are performed not on the full matrix element but in the leading 
pole approximation of soft and collinear regions in 4 + e dimensions. The soft and 
collinear poles now appear as poles in e and in addition the soft part would depend 
logarithmically on the soft cut-off 6s while the collinear part would depend on the 
both Sg and 6c- All positive powers of the small cut-off parameters are set to be zero. 
The phase space degrees of freedom that remain, correspond to a 2-body process 
and can now be combined with the virtual diagram. 

The virtual corrections to the annihilation process to order 0{as) can be obtained 
by considering the gluonic virtual corrections to the vertex and wave function renor- 
malisation for the process gg — )■ Z and then attaching the KK modes to all possible 
legs and vertex as allowed by the Feynman rules [21 II]- This would generate 27, one 
loop 2 — 2 diagrams which on multiplication with the q q ^ Z Gkk at LO would 
give all the virtual contributions for the annihilation process. Performing the loop 
integrals in 4 + e dimensions would lead to poles in e in the soft and collinear regions 
and combining this with the real diagrams would lead to the cancellation of the soft 
singularities. The remaining collinear singular terms which appear as poles in e are 
systematically removed by collinear counter terms in the MS factorisation scheme, 
at an arbitrary factorisation scale fip- The resultant expression is now finite but 
depends on the cut-off parameters 6s and 6c logarithmically. Combining this 2 — )■ 2 
part with the finite 2 — )■ 3 part and performing the phase space integration for the 
various distributions, it is expected that the results be independent of the choice of 
the slicing parameters 6s and 6c- 

The q{q) g — > q{q) Z Gkk process begins at 0{as) so does not get any virtual 
corrections to this order and can be obtained from the 2 — i- 3 annihilation diagrams 
by using the crossing symmetry. The analytical results of LO 2 — )■ 2, finite part of 
NLO virtual corrections and full 2 — )■ 3 matrix elements that go into our numerical 
code will be presented in the longer version [2T]. 

3 Numerical results 

In this section we present the numerical results for the associated production of the 
Z-boson with the KK modes at the LHC {y/S = 14 TeV) to NLO in QCD. The 
mass of the Z-boson and the weak mixing angle are taken to be mz = 91.1876 and 
sin^^^ = 0.2312 respectively. The fine structure constant is taken to be a = 1/128. 
Through out our computation we have used CTEQ6 L/M parton density sets [22] 
and Hf = 5 light quark flavors. The two loop running strong coupling constant in 
the MS scheme has been used with the corresponding Aqcd = 0.226 GeV. Unless 
mentioned other wise, both the renormalization and the factorization scales are set 
to iiR = hf = Pt, the transverse momentum of the Z-boson. Further, the following 
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cuts have been implemented in our numerical code: 

p|, p^''' > p^'"" and \y^\ < 2.5 (5) 

where is the rapidity of the Z-boson, p™" = 400 GeV and the missing transverse 
momentum p^**** is given by 

p™- = p| (pG) for pP;t < 20 ( > 20) GeV (6) 

At LO, the missing transverse momentum p™*** is the same as pf . The additional 
jet at NLO can be soft or hard; For the soft jet the p™*** is the same as pf while for 
hard jet pjp^'^^ is p^? of the KK modes (Eq. ([6])). In addition to the above, we have 
put a cut on the pseudo rapidity of the jet \ri^'^^\ < 2.5 for p^^* > 20 GeV. 
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Figure 1: Dependence of the order contribution to the transverse momentum 
distribution of the Z-boson at the LHC, on the slicing parameter 6s (top) with 
= Ss/100 and for = 3 TeV and 6 = 4. Below the variation the sum of 2-body 
and 3-body contributions is contrasted against the value at Sg = 10~^. 

We check the stability of our results against the variation of the slicing parameters 
6s and 6c introduced in the slicing method. In Fig. [H we present the dependency 
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of both the 2-body and 3-body pieces of the 0{as) contribution on 6s, keeping the 
ratio Sg/Sc = 100 fixed. These results are obtained for a specific choice of the ADD 
model parameters M^, = 3 TeV and 6 = 4. It can be seen from the figure that 
the sum of these two pieces is positive and is fairly stable for a wide range of 6s- 
The positive sum implies that the QCD corrections do enhance the leading order 
predictions. In the rest of our work, we choose 6s = 10~^ and 6c = 10~^. 
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Figure 2: Transverse momentum distribution of the Z-boson (left) and the corre- 
sponding K-factors (right) for the associated production of the Z-boson and the KK 
mode at the LHC, for = 3 TeV and 5 = 2,4, 6. 

In Fig. [2], we present the transverse momentum distribution of the Z-boson at 
the LHC to NLO in QCD and its dependence on the number of extra dimensions 6 
for = 3 TeV. In each of the distributions corresponding to 5 = 2, 4, 6 the QCD 
effects are found to have increased the leading order predictions considerably. In 
this transverse momentum distribution of the Z-boson, the i^'-factor, defined as the 
ratio of NLO cross sections to the LO ones, is found to increase with and vary 
from 1.1 to 1.46 depending on the number of extra dimensions 6. In a similar way, 
the missing transverse momentum distribution is shown in the left panel of Fig. |3l 

The ADD model which is an extension of the SM to address the hierarchy prob- 
lem is an effective theory and the UV completion of the TeV scale gravity has to 
be quantified. For the real KK mode production process, the kinematical constraint 
discussed earlier, provides a natural UV cutoff on the integration over the n-sphere, 
but the hard scattering scales involved at the LHC energies can be close to the fun- 
damental scale Mf). To study the sensitivity of our results close to the UV region [3], 
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Missing transverse momentum distribution 
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Figure 3: Missing transverse momentum distribution (left) at tlie LHC for Mo = 
3 and for 6 = 2 and 4. Tlie truncated as well as the un-truncated transverse 
momentum distributions of the Z-boson (right) at the LHC for Md = 5 and for 
(5 = 4 for both LO and NLO. 



we compare at LO and NLO the pt distribution of Z-boson wherein the invariant 
mass of the KK mode and Z-boson Qzg is (a) computed only when Qzg < Mo 
(truncation) and (b) computed for all possible values of Qzg (un-truncation) . In 
the right panel of Fig. [3l we compare the results for the ADD model parameters, 
Mo = 5 TeV and (5 = 4 at LO and NLO for the pt distribution. As compared to 
the un-truncated distribution, the percentage difference is tabulated below for LO 
and NLO for a few values of . 
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The difference between the un-truncated and truncated results become larger with 
(a) increase in p^, (b) increase in the number of extra dimensions 6 and (c) decrease 
in the fundamental scale Mo- Further it is seen that the NLO QCD corrections do 
decrease the difference between the un-truncated and truncated results. 

Finally, we study the dependence of both LO and NLO cross sections on the 
factorization scale /i_F by varying it from 0.2 p§, to 2.0 . One of the motivations 
for the computation of the QCD corrections is to minimise the scale uncertainties 
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Figure 4: Scale uncertainty in the transverse momentum distribution of the Z- 
boson at the LHC (left), for a variation of the factorization scale [ip vn the range 
[0.2, 2.0] pf and for the choice of Md = 3 TeV and 6 = 4. The transverse momentum 
distribution of the Z-boson for = 7 TeV at the LHC (right). 



by computing the cross sections to higher orders in the perturbation theory. As 
expected, the scale uncertainties in the leading order predictions are considerably 
decreased after incorporating the one-loop QCD corrections to the associated pro- 
duction of the Z-boson and the KK modes at the LHC. The results are shown in the 
Fig. m (left) for the case of transverse momentum distribution of the Z-boson, for 
Md = 3 TeV and 5 = 4. In Fig. H] (right) we have also plotted the pt distribution 
for the current LHC energies of ^ = 7 TeV. The K-factor ranges from 1.05 to 1.13 
for 500 < < 1000 GeV. The contribution of q{q) g subprocess for the = 7 
TeV is much smaller than the contribution at -s/S* = 14 TeV and that accounts for 
the much lower K-factor. 

It is interesting to note here that a similar analysis goes through for the associated 
production of and the KK modes in the large extra dimensional model at the 
LHC. The difference lies both in the couplings of the quarks to the weak bosons and 
in the respective parton fluxes to be convoluted with the partonic cross sections. In 
Fig. [5] we have plotted the W~ distribution for = 14 TeV at LO and NLO. For 
the distribution the K-factor is in the range 1.25 - 1.37 for 500 < p^ < 1000 
GeV. A complete analysis of the associated production of in association with 
KK modes of the ADD model will be presented in the longer version j21] . 
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Figure 5: The pr distribution of W in association with KK modes at the LHC 
with VS = U TeV. 

4 Conclusion 

To conclude, we have computed the NLO QCD corrections to the associated produc- 
tion of the vector boson at the LHC, using semi-analytical two cutoff phase space 
slicing method. We have presented results for = 14 and 7 TeV. Our results are 
checked for the stability against the variation of the slicing parameters 6s and 6c- 
We have studied the truncated as well as the un-truncated transverse momentum 
distributions of the Z-boson, together with the missing transverse momentum distri- 
bution and their dependence on the number of extra dimensions 6. The NLO QCD 
corrections are found to have not only enhanced the LO cross sections considerably, 
with the K-f actors ranging from 1.1 to 1.46 depending on the 6 = 2, 4, 6 for Md = 3 
TeV, but also decreased their factorization scale uncertainties significantly. 
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